Effective density dependent pairing forces in the T=l and T=0 channels. 
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Effective density dependent pairing forces of zero range are adjusted on gap values in T = 0, 1 
channels calculated with the Paris force in symmetric nuclear matter. General discussions on the 
pairing force are presented. In conjunction with the effective k-mass the nuclear pairing force seems 
to need very little renormalization in the T = 1 channel. The situation in the T = channel is also 
discussed. 
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I. INTRODUCTION 

The novel availabihty of exotic nuclei has spurred 
an immense revival of nuclear structure investigations 
[0. Indeed nuclei close to the neutron or proton drip 
lines may exhibit very unusual features such as pro- 
nounced neutron or proton skins and neutron ha- 
los Q. Among many very interesting questions, nuclear 
pairing has again become on the forefront of theoretical 
interest. Indeed the existence of neutron halos is due to 
the pairing force and in heavier proton rich N Z 
nuclei the proton-neutron pairing may play an important 
role 1^. In this work we therefore want to address some 
problems of neutron-neutron and proton-neutron pairing. 
This concerns for instance considerations of the effective 
pairing interactions. However, we also will discuss some 
other aspects of more general character. We mostly will 
study the infinite matter case. 



II. GENERALITIES ON THE NUCLEAR 
PAIRING FORCES 



It is a well established fact that, aside from the excep- 
tion of magicity, nuclei are superfluid. There can be nn 
as well as -p-p pairing whereas "pn pairing is less frequent. 



One of the main questions we will treat here is the ef- 
fective pairing force. We will do this in the framework 
of homogeneous nuclear matter at various densities. The 
limit to finite nuclei can be established through the Lo- 
cal Density Approximation (LDA) which seems to work 
very well also for the nuclear pairing problem j7j. Quite 
generally the equation for the gap A in nuclear matter 
can be written as 



A, 



(1) 



where Wpk is the (effective) pairing force, the are the 
Brueckner-Hartree-Fock single particle energies and is 
the Fermi energy. The summation goes over momentum 
states. In (|l|) we did not specify whether we consider the 
T = 1 or T = channels. 

The first aspect we want to discuss is what kind of force 
fpfe shall be used in Eq.(0) from a microscopic point of 
view. The answer to this question is in principle very 
well known since the early days of superconductivity and 
superfiuidity. Since the gap equation can be derived from 
the Bethe-Salpeter equation for the two particle many- 
body Green's function ||,^ , the pairing force v^k is built 
out of the sum of all particle-particle irreducible Feyn- 
man graphs . To lowest order in the bare interaction 
it is given by Fig. 1. 





FIG. 1. Schematic representation of the pairing force v^h to lowest order in the bare interaction. 



In Fig. 1 the dot stands for the bare vertex. The 
second term represents a ph screening correction to the 
bare force. The very important point we want to make 
here is that in no way some type of Bethe-Goldstone or 



Brueckner G-matrix can be used in the gap equation. 
Since the gap equation is already a kind of in medium 
two-body Schroedinger equation (see e.g. Ref. [ p^ , pd] |) 
one cannot use a G-matrix which in itself is a solution 
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of the in medium two-body problem. Otherwise there is 
severe double counting. Since the G-matrix essentially 
softens the short range repulsion one expects that pair- 
ing becomes enhanced if used in the gap equation. In the 
pairing problem everything depends exponentially on the 
system parameters 12 1 and this effect can then be quite 
large. A demonstration is given in Fig. 2, where the 
nn gap is calculated once with the bare Paris force 
and once with the corresponding G-matrix jl^ . One sees 
that the use of the G-matrix enhances the gap value by 
practically a factor of two. 
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FIG. 2. Pairing gap Af in neutron matter as a function 
of the Fermi momentum kp with the Paris force and with the 
corresponding G-matrix. 

Sometimes Eq.(^ is divided into a low momentum and 
a high momentum space and the high momentum space 
is eliminated in renormalizing consistently the bare in- 
teraction in the low momentum space |I4|. This type 
of procedure is, of course, perfectly allowed, since it is 
only a mathematical trick for solving Eq.(^. Unfortu- 
nately in nuclear physics it is a quite widespread habit 
since decades (see for example and the critics given 
in ) to use some kind of G-matrix in Eq. ([^) as for ex- 
ample Skyrme forces which are to be considered as a phe- 
nomenological representation of a microscopic G-matrix. 
One will object that one of the most successful nuclear 
nn pairing forces, namely the Gogny force is also to 
be considered as a G-matrix. Things are, however, more 
subtle there as we now will explain. The first observa- 
tion is that the Gogny force in the ^ 5*0 channel is of finite 
range but density independent. Second, one finds when 
solving the gap equation with the Gogny force in nuclear 
matter that it gives results which are very close to the 
ones obtained with the Paris force or any other realistic 
bare nucleon-nucleon force. This is demonstrated in Fig. 
3 where we compare results of the gap from the DI, DIS, 
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FIG. 3. Pairing gap Af in the ^So channel in symmetric 
nuclear matter calculated with the Gogny forces Dl and DIS, 
compared with the Paris force results, from Refs. 



We see that DIS is still much closer to Paris than Dl. 
Indeed DIS has been readjusted ||r^ to give in first place 
a lower surface tension than Dl but at the same time 
to give a smaller even-odd staggering so that it becomes 
in closer agreement with experiment. It is very surpris- 
ing that this readjustment brought DIS so close to the 
bare Paris force. So in the S = T = 1 channel the 
Gogny force acts like a realistic bare force at least in 
what concerns energies up to the Fermi energy. This 
conclusion was also found in Q and is further confirmed 
by the fact that the scattering length corresponding to 
DIS, aois — 12.16 fm, is very large and of the same 
order of magnitude as the experimental value a ~ 18.5 
fm. 

The reason why the Gogny force acts like a free force 
in the nn pairing channel in spite of the fact that it has 
been adjusted to the G-matrix from the Sprung- Tourreil 
force [|l8[ can only be guessed: probably for this force in 
that channel the Pauli blocking is so efficient that in the 
G-matrix equation, G — v + v^G , the second term on 
the r.h.s. is suppressed. On the other hand the question 
remains why experiment apparently demands a pairing 
force very close to the bare one. This is true at least 
in the T — 1 channel. For the T = channel much 
less investigations have been performed and it is unclear 
whether a bare force can be used as well. One reason 
which can be advanced to explain the validity of the bare 
force is a possible cancellation between screening effects 
and effective mass enhancement. Graphically these two 
possibly opposing effects are shown to lowest order in the 
interaction in Fig. 4. 
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FIG. 4. Schematic representation of the effective mass enhancement and of the screening effect to lowest order in the 
interaction. 



In this respect it should be mentioned that the Hartree- 
Fock-Bogohubov (HFB) calculations with the Gogny 
force are performed with the so called fc— mass m* < m. 
However one knows that the corresponding level density 
close to the Fermi energy is much too small. Includ- 
ing i^— mass corrections such as the one shown in Fig. 
4 brings the effective mass at the Fermi level back to 
the bare mass or even overshoots it. For consistency the 
screening of the bare force also shown in Fig. 4 must be 
included. Larger effective masses enhance pairing while 
screening probably weakens it so that the net effect could 
be the bare force. To investigate such effects, extreme 
care must be taken that both contributions of Fig. 4 
are treated on the same footing. Since, as already men- 
tioned, pairing depends exponentially on the system pa- 
rameters, the slightest imbalance (for example in treating 
both graphs of Fig. 4 in slightly different approxima- 
tions) may cause strong erroneous results. One way to 
treat things consistently could be to use the Gorkov equa- 
tions p9| and develop the normal and abnormal parts of 
the mass operator matrix to second order Born approxi- 
mation and solve the corresponding gap equation numer- 
ically. In medium effects similar to the ones shown in Fig. 
4 have been included in the past to the pairing problem 
in one way or the other ||2^. Practically all calculations 
resulted in an important reduction of A = A (kp) com- 
pared to the values shown in Fig. 3. It can be deduced 
from the study in that a reduction of pairing in infinite 
matter obtained with the Gogny force in a global way, i.e. 
for all values < fci;- < 1.4 fm~^, inevitably leads also to 
a reduction of pairing in finite nuclei of the same propor- 
tions (this fact can be understood via the local density 
approximation which as mentioned already, on average, 
yields comparable results to quantal calculations |]7|,pl[). 
It therefore can be concluded that e.g. a reduction of the 
A— values in Fig. 3 by a factor of two (a scenario often 
encountered in the calculations of references mentioned 
above) will fail to reproduce experimental gap values of 
nuclei when the underlying theory is applied to finite nu- 
clei. 

Concluding these general considerations we want to say 
that in absence of any necessity stemming from experi- 
mental facts it is probably safe to treat nuclear pairing 
in conventional mean field theory with the bare nucleon- 
nucleon potential as this is indicated from the micro- 
scopic theory and as apparently is needed to reproduce 
experimental facts in the T = I channel. Using this phi- 
losophy one arrives naturally for T — np pairing at 
much stronger gap values since the N — N force is 
strongest in this channel. We will give some more details 



about this in the next section and also discuss how the 
bare interaction in the gap equation can be replaced by 
an equivalent density dependent zero range force such as 
they have become quite popular recently in the nuclear 
structure problem. 



III. EFFECTIVE DENSITY DEPENDENT ZERO 
RANGE PAIRING FORCES 



In the last section we gave arguments that, at least as 
a first guess, it is indicated to use as the pairing force the 
bare nucleon-nucleon potential. We here want to develop 
arguments that this strategy is not necessarily orthog- 
onal to the popular employment of density dependent 
zero range forces with a cutoff. Such arguments have 
first been developed by Bertsch and Esbensen Q and we 
here want to refine these arguments on the one hand and 
on the other side extend them also to T — np pairing. 
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FIG. 5. Vkpkp vs. kp in the 5 = T = 1 channel for the 
Gogny Dl force. 

A qualitative argument why a density independent fi- 
nite range force in the gap equation (Eq.(|lj)) can be re- 
placed by a density dependent zero range one with a cut- 
off, goes as follows. For s— wave pairing only the an- 
gle averaged matrix element Vpk enters the gap equation 
Ap = '^i^VpkKk, where Kk = Ak/2Ek is the abnormal 
density and 
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Eh — 



(2) 



the quasiparticle energy. The former is very much peaked 
at k — kp with a peak width of the order A = Afe^ . Since 
anyway in pairing problems only the gap values at fc ~ kp 
matters, we see that for A^.^ only the value of the ma- 
trix element Vkj,kF plays a significant role. In the Gogny 
force, this matrix element as a function of hp is shown in 
Fig. 5. Since a 5— force is a constant in fc— space, one has 
to weight the 5— force with a kp, i.e., a density dependent 
factor similar to Vkpkp in order to recover the essential 
pairing features of the original finite range force. The 
only thing we have to add is a cutoff value, otherwise the 
gap equation would diverge. Bertsch and Esbensen 
therefore proposed the following density dependent zero 
range force 



w(i"i,r2) 



vo 



l-rj 



PO 



Sin-r^) (3) 



where wo,?y, a are adjustable parameters and po is the 
saturation density. In the gap equation (|l|), Eq.(H) must 
be supplemented with a cutoff value ec which thus con- 
stitutes a fourth parameter. However, at zero density 
the cutoff and vq must be chosen such that the scatter- 
ing length a is reproduced. For Eq. (^) one obtains the 
relation 



Vo 2tt^ 

m 



+ kc 



(4) 



where — ec- The neutron-neutron scattering length 
is very large (a — 18.5 fm) and we take in Eq.(U) the 
limit a oo that leads to the relation for vq also used in 
Q . One finally remains with three adjustable parameters 
(77, a, ec)and the gap equation reads 



1 = 

_Vq_ 

,r2 



(5) 



Po 



2r 



de 



{e-epy + A'^ 



With respect to Q we considered also a density depen- 
dent effective mass. Since finite nuclei calculations are 
performed with such an effective mass one must account 
for it when adjusting a 5— force which later shall be used 
in BCS or HFB calculations. For the effective mass we 
take the one corresponding to the Gogny force. 
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FIG. 6. T — 1 pairing gap in nuclear matter. The dots 
are the results of a Hartree-Fock calculation using the Gogny 
force. The continuous (a) and dashed (b) curves are the re- 
sults obtained with the effective pairinginteraction in Eq.(^) 
with (a) effective mass m*/m as in Eq.M and (b) m* /m — 1 
(see text). The dotted line corresponds to the pairing force 
in Ref. @ (see text). 

In Fig. 6 we show the fit to A (kp) in the isovec- 
tor channel obtained from Eq. (||) with ec = 60 MeV, 
77 — 0.45, a — 0.47 . Also shown is the fit corresponding 
to the bare mass (i.e., m* /m ~ 1) with ec = 60 MeV, 
•q — 0.70, a = 0.45 , as in Ref. Q. In both cases, the cor- 
responding Vq value is vq =481 MeV fm"^. We see that 
the fits are good for values of kp up to the saturation 
value kp = 1.35 fm^^. A density dependent (5— force has 
also been used for T = 1 pairing in finite nuclei in the 
context of the HFB pal and in the context of relativistic 
Hartree-Bogoliubov [^3| . The strength used there is how- 
ever larger. If we use the pairing force in Ref. ||2^] with 
Vo = 700 MeV fm^, we get the dotted line curve shown 
in Fig. 6 that corresponds to the following parameters in 
our notation: vo = 1400 MeV fm^; ec = 7 MeV; ?7 = 1 
MeV and a = l MeV. 

For finite nuclei, the force 
proximation 



can be used in BCS ap- 



(m \v\ kk) 



k,£k<<^C 



Afc 
2Ek 



(7) 



m* {p)\ ^ , m kp ^ 



1 + + 2 {B, - H,) - AM, 



Ple~^^ 



cosh(xc) sinh (a;c 



(6) 



or in the HFB approach where the gap equation has the 
form (|^) in the canonical basis. We want to point out 
that the cutoff has to be counted relative to the bottom 
of the single-particle well and not from its edge. 



with Xc = kpfil/2, and the coefficients WpBc, H,, M,, fJ,c 
corresponding to the Gogny force Dl [p|,p^. 
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IV. PROTON-NEUTRON PAIRING IN THE T= 
CHANNEL 



In this section we want to extend our considerations 
to n — p pairing in the T — 0, i.e. in the deuteron chan- 
nel. As we suggested earher, as a first guess one should 
investigate the gap equation with the bare force. The 
gap equation in homogeneous symmetric nuclear matter 
has recently been solved for the T = channel using 
the bare Paris force with single particle energies obtained 
in Brueckner-Hartree-Fock approximation. Since in the 
deuteron channel (T = 0, 5 = 1, L = 0, 2) we have a mix- 
ture of s— and d— waves involving the tensor force, the 
net outcome is more attraction leading to the deuteron 
bound state in free space. This increased attraction then 
takes over to the gap equation (which in the zero den- 
sity limit turns into the Schroedinger equation for the 
deuteron, see [^,|^) and, not unexpectedly (remember 
the exponential dependence), the gap values in the T = 
channel as a function of kp are much stronger reaching 
values more than a factor of two larger than in the T = 1 
channel. This is shown in Fig. 7 (Ref. pT[|). 
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FIG. 7. Pairing gap versus Fermi momentum for sym- 
metric nuclear matter in the T — channel from the Paris 
potential. 

The use of the bare force in the T = channel may, 
however, be more questionable than in the T = 1 chan- 
nel. This stems from the implication of the d— wave, i.e. 
the tensor force. The latter seems to be more affected 
by medium effects than the s— wave part and therefore 
certainly great care must be employed in this channel. In 
particular, it has been shown in ||25|| that higher shell ad- 
mixtures make the tensor force appear weaker in the va- 
lence space. Again the possible balance of the two graphs 
of Fig. 4 should thoroughly be investigated with respect 
to s— and d— wave contributions. We do not exclude 



the possibility that the tensor force is largely screened 
in the medium and thus the enhancement of the T = 
gap values may be brought back closer to values to which 
we are used in the T = 1 case. However, without hav- 
ing detailed investigations at hand, we here stick to our 
working hypothesis and base our considerations on the 
bare force scenario. In this sense it may be interesting 
to also adjust, like we have done it for the T = 1 case, 
a density dependent ^— force to the calculation with the 
Paris force shown in Fig. 7. In principle, in this case, the 
parameter vq should be chosen such that the deuteron 
binding energy is reproduced in free space. We, however, 
found that with this condition the cutoff parameter must 
be chosen very large rendering this force not very prac- 
ticable in actual calculations. We therefore adopted the 
strategy to also vary within very narrow limits the pa- 
rameter vo what may slightly degrade the gap values at 
very low densities but significantly improves them at the 
higher densities. In Fig. 8 we show such an adjustment 
using various cutoffs. The value of vq used for the fits in 
Fig. 8 is 



-1.05 



m kc 



These fits should be useful for finite nuclei calculations. 



• Paris 

X e_,=60 MeV; a=0.57; ri=0.54; bare m 

- E^=80 MeV; a=1 .00; ri=0.25; m* 

- e^=60 MeV; a=0.75; ri=0.20; m* 

- e^=50 MeV; a=0.70;ri=0.15; m* 

- e^=40 MeV; a=0.70; ri=0.10; m* 

- E^=30 MeV; a=0.70; ri=-0.10; m* 




k,(fm-) 

FIG. 8. r = pairing gap in nuclear matter. The dots 
are the results obtained from the Paris potential. The various 
curves correspond to fits with Eq.(|^, using difi^erent param- 
eters. 
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V. CONCLUDING REMARKS 

In this paper we critically reviewed the use of effec- 
tive nuclear pairing forces. We argued that a Bethe- 
Goldstone or Brueckner G-matrix must not be used in 
the gap equation. As a first guess, not knowing any- 
thing better, the free nucleon-nucleon force may be tried 
in the gap equation. At least in the traditional T = 1 
channel this prescription seems to work remarkably well, 
since the best phenomenological force, namely the Gogny 
force, acts very nearly like a free force in the T — 1 pair- 
ing channel. We then advocated that the same strategy 
should be adopted in the T = channel. We pointed out 
that the situation may, however, be slightly more subtle 
there because it is the action of the tensor force which 
makes the T = channel more attractive than the T = 1 
one. The tensor part of the nuclear interaction is, how- 
ever, a very delicate subject and it may well be that it 
is more affected by screening than the rest of the force. 
In the second part of the work we demonstrated that the 
use of density dependent zero range forces in the pair- 
ing channel may not be orthogonal to the use of finite 
range density independent forces. Following Bertsch and 
Esbensen |^ , we give parametrizations of density depen- 
dent ^— forces which reproduce the gap values in both 
T = and T — 1 channels very well over the whole 
range of relevant nuclear matter densities. Such forces, 
augmented by a cutoff, should then also be useful for 
calculations in finite nuclei. 
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